On superstatistical multiplicative-noise processes 
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In this article we analyse the long-term probability density function of non-stationary dynamical 
processes with time varying multiplicative noise exponents which are enclosed inwards the Feller class 
of processes. The update in the value of the exponent occurs in the same conditions as presented 
by Beck and Cohen for superstatistics. Moreover, we are able to provide a dynamical scenario for 
the emergence of a generalisation of the Weibull distribution previously introduced. 



I. INTRODUCTION 



The description within a physical context of driven 
non-equilibrium complex systems has frequently been 
made by considering that their dynamical behaviour is 
characterised by spatial-temporal fluctuations of some 
parameter, (3. Usually, this parameter has been con- 
sidered to be the (inverse) temperature, the dissipation 
of energy in turbulent flows, the amplitude of Gaussian 
white noise, the local mean-reverting value or the lo- 
cal variance. As an example, we mention the standard 
case of a Brownian particle diffusing along an inhomo- 
geneous medium in which temperature (hence diffusion 
"constant" ) fluctuates in both space and time. In this 
approach, as it can be understood, there are two impor- 
tant time scales: the scale in which the dynamics is able 
to reach a stationary state (assuming a fixed value for 
parameter 0), and the scale at which the fluctuating pa- 
rameter evolves. A particular case to consider is when 
these two time scales are clearly separated, specifically, 
when the time needed for the system to reach stationarity 
(considering a predetermined 0) is much smaller than the 
scale at which that parameter changes. In the long-term, 
the non-equilibrium system is described by the superposi- 
tion of different local dynamics at different time intervals 
that was coined by Beck and Cohen as superstatistics 
or "statistics of statistics" P, @|. Frequently, systems 
that are characterised as "superstatistical" exhibit non- 
Gaussian distributions with kurtosis excess, or distribu- 
tions with non-exponential decay. In addition, supersta- 
tistical systems present a parameter, /3, that fluctuates 
on a large scale, T, and follows a time-independent dis- 
tribution, p{0). The superstatistical framework has suc- 
cessfully been applied on a widespread of problems like: 
interactions between hadrons from cosmic rays fluid 
turbulence JEjE B, g ranular material 0] , electronics § , 
economics [|[ flfl TlilfT^| , among many others [HI . Fur- 
thermore, it has been regarded as a possible foundation 
for non-extensive statistical mechanics [3| based on Tsal- 
lis entropy [l4j as we show later on. In this manuscript 



we introduce a different analysis of differential stochas- 
tic dynamics in which the exponent of a Feller process is 
assumed as the superstatistical parameter. The main ad- 
vantage of this proposal is that it permits the evolution of 
the functional form of the second-order Kramers-Moyal 
moment in opposition to previous presentations. 



II. SUPERSTATISTICS 

Consider an inhomogeneous system composed by a 
large set of cells that have different values of some pa- 
rameter f3 as we have referred to here above. Within each 
cell, local equilibrium is reached very promptly. The pa- 
rameter (3 is taken as constant throughout a period of 
time T after which it changes into a new value. This up- 
date occurs always in accordance with a distribution p((3) 
for the parameter. Taking into consideration that each 
cell is in local equilibrium, thus presenting a Boltzmann 
factor, e~ l3E ' , x the long-term stationary distribution 2 of 
the non-equilibrium system is obtained from a weighted 

average of local Boltzmann factors, with (3 = (3 ^/3^ , 



P(E) = pip) p(E) 



e -0B 
W) 



d0, 



(1) 



where p (E) represents the density of states, and Z (0) 
the normalisation constant. Going back to the example 
of a Brownian particle moving across a medium treated 
in Ref. [|[, we get that its velocity, v, is obtained from 
the local Langevin equation, 



d v = —7 v dt + a dWt 



(2) 



Seeing that the medium is inhomogeneous, either 7 [151 ] 
or a vary from cell to cell on a large time scale T 3 . 
Therefore, local Boltzmann factor has, 



= 



2-/ 



(3) 
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E is the effective energy in each cell. 

2 The observation time t^>T. 

3 If 7 is the random parameter then, /3 = 7, else if o is the random 
parameter then, j3 = cr. 
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which is random. From Eq. ((3|) , the parameter (3 can also 
fluctuate for the case of a particle with varying mass [l6[ . 

Within time scale T, and according to Eq. ((2]), the 
local stationary distribution of velocities is a Gaussian 
conditioned to value /3, 



p' (v\P) = 
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(4) 



If the system is able to reach some local equilibrium 
before an update of $ takes place, i.e., T ;» 7 _1 = r, 
then, we can determine the marginal velocities probabil- 
ity distribution of the long-term behaviour of the Brow- 
nian particle, 



P(v) = / P (/3)p'(v\P) d(3. 
Jo 



(5) 



Hence, it is straightforward to verify that the form of 
P (v) depends explicitly on the functional form of p (/3) . 
Specifically, it was verified in Ref. [3j that, when p(/3) is 
the x 2 — distribution with n degrees of freedom, Eq. ([5]) 
yields, 

P{v) = ^[l + {l- q )^] 1/{l ~ q \ (6) 

where q = 1 + ^xjj Z is the normalisation factor, and (3q 
the average inverse temperature (see Ref. [l| for details). 
Such a distribution P (v) maximises Tsallis entropy |14| , 



1-J\p(x)] q dx 



{q G K) 



(7) 



This fact has turned out superstatistics into the first dy- 
namical scenario for the emergence of non-extensive sta- 
tistical mechanics [13] • 



III. THE MODEL 



whose solution / (v,t) relaxes exponentially with a char- 
acteristic time, r, into the stationary solution, 



P (v) = ^ exp 
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cu 2 (l-2a 
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2(l-2a) 



(10) 



i.e., a Weibull-like distribution, W(v). Z is the normal- 
isation constant, J p(v) dv, 



1 -4a 
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u 2 {l-2a) 
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(11) 



For a = 0, Eq. (JSJ) becomes the standard Langevin equa- 
tion, and p (v) the Gaussian distribution, 



g(v) 



i 



■ exp 
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with Z 



7TU^ 
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After the transient, f (v,t) w The mean value 

of w, £ = f vp (u) du, is equal to zero as well as all odd 
moments of p(v). Regarding the second-order moment, 
v 2 = J v 2 p (v) dv, we have got 
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if a = 
if a ^ 



If we consider v as the velocity of a particle of unitary 
mass which does a ID random walk, using equipartition 
theorem we are able to determine the inverse tempera- 
ture, (5 = (fcT)" 1 5 , yielding 



(3 = v 2 



(14) 



Evaluating the kurtosis 



Consider the following one-dimensional stochastic dif- 
ferential equation, 



dv 



--fvdt + u [v 2 ] a dWf 



{v ^ if a < 0) . 



(8) 

where Wt is a regular Wiener stochastic process, i.e., 
(dW t ) = 0, and {dW t dW t >) = dt8{t-t') 4 . Stochas- 
tic equation ([8]) belongs to the Feller class of (multi- 
plicative noise) processes [HI with 7 > 0, a < \ for a 
(time-dependent) normalisable probability density func- 
tion (PDF) f(v,t). The associated Fokker-Plank Equa- 
tion of Eq. ([8]) is 



df(v,t) d 



dt 



dv 



[jvf(v,t)} 



28^ 



f(v,t) 



(9) 



v 4 p (v) dv, 



we have obtained 
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As it is visible from Fig. [T] distribution (fT0| is platykurtic 
for a < 0, and leptokurtic for a > 0. 

Moving on, we shall now consider that exponent a in- 
stead of constant, varies according to superstatistical re- 
quirements. In other words, let us consider an inhomo- 
geneous system which is composed by a large set of cells 
that have different values of a. Within each cell, the 



4 The main advantage of writing [f 2 ] Q instead of 
a' = 2a is that of analyticity for all v when a > 0. 



with 



In this expression T is the tempetature. 
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FIG. 1: Kurtosis, k, vs. a. Upper panel: Non-positive 
values of a which lead into platykurtic distributions. In the 
limit a — > — oo, k = 1. Lower panel: Non-negative values of 
a leading into leptokurtic distributions which has a = | as 
upper bound. The dashed line corresponds to the kurtosis of 
a Gaussian k, — 3. 



value of a is updated at every elapsed time interval £, in 
agreement with a certain PDF p (a) . The update scale, 
£, is much greater than relaxation time scale r. Notic- 
ing that each cell is in local equilibrium, the long-term 
stationary distribution of the non-equilibrium system is 
obtained performing the integral, 



P ( v ) — / p(ot) p{v) da. 



(17) 



As possible applications of such a model we name: de- 
scription of velocities in granular material (particularly 
see figures of Ref . [19j ) and unconventional turbulent flu- 
ids, or even the dynamics of financial observables. Ex- 
plicitly, systems which undergo through different phases 
during their time evolution or situations where different 
stages are measured when observations are made at the 
same point of space. In a thermodynamic context, the 
fluctuations in a correspond to fluctuations in temper- 
ature, but obtained through a completely different way 
from the proposal presented in Ref. Q . Specifically, fluc- 
tuations in a induce a modification of the functional 
form of the 2 nd order Kramers-Moyal coefficient, while 
in Ref. 0] its functional form is always preserved. 



A. Some examples 

1. Dichotomous case 

This case represents the simplest form to introduce 
fluctuations in a, and for which a full analytical treat- 
ment is possible. Its probability density function is sim- 
ply, 



P (") = t; s ( a - a o) + 7;$ (« ~ oil) ■ 



(ao^ax). (18) 



Amongst all endless possibilities for cvo and a\, let us 
firstly consider cases for which one of the exponents is 
equal to zero. The long-term distribution is thus given 
by 



P{v) = \Q{x) + \w a {x). 



(19) 



For small values of \v\, and a < 0, P (v) approaches the 

-2a 



limit -yj - 
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In Fig. [2] we exhibit the resulting probability density 
function, P(v), for ap = 0, and a% = — 1, ax = — |, and 



a = I. The value v^ f , 



u ef 



v 2 P (v) dv, 



(20) 



is presented in Fig. [3] for several values of a\ 
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FIG. 2: Probability density function P (v) vs. v for the case 
of superstatistical processes with a dichotomous distribution 
(fT8)l where 7 = u = 1. Upper panel: In both cases ao = 0, 
and cti — — 1 for the full line, and ai = — | for the dashed 
line. Lower panel: P (v) vs. v for a dichotomous case with 
ao = and a = | (7 = w = 1) in a log-log scale. For this case, 
the variance v\j — 0.2860 . . . and kurtosis n — 11.953 .... 
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FIG. 4: Numerically obtained probability density function 
P (v) vs. v for the case of superstatistical processes with uni- 
form distribution, Eq. (|21|l where 7 = u — 1. The black line 
corresponds to Qo = 0, and ai = — |, with = 0.833..., 
and kurtosis, k = 1.755 . . .. The grey line corresponds to 



a uniform distribution with 



no 



ai = 



yielding 



"ef 



1.169. 



and kurtosis, 



1.179. 



FIG. 3: Standard deviation, irj., vs. dynamical exponent 
ot\ for superstatistical process with a dichotomous PDF fj 18p 
where ao = 0, and 7 = w = 1. For on = 0, «^ = |, 
and for ai = —00 the value of v^. tends to | (upper panel). 
For positive values of ai, v\* is a monotonically decreasing 
function of ai with = | for m = \. 



3. The x -distribution 

Another distribution that appears in various phenom- 
ena is the ^-distribution, 



p(a) 











r 1 






j exp 


--H] 






L a i 



(22) 



2. Uniform distribution case 

If, for the previous example the analytical form of P (v) 
is easily obtained, that does not happen in situations 
in which a evolves according to a uniform distribution 
between ao e a±, 

— - — if ai < a < an 

an— ai L — — u 

/'(Ml { . (21) 

otherwise 

However, it is possible to evaluate numerically the form 
of the long-term distribution as we present for two par- 
ticular cases in Fig. |H Computing the kurtosis, n, for 
both cases we have verified that the two examples are 
platykurtic. 

Another possibility is to consider just non-negative val- 
ues for a. This situation is more likely to be experimen- 
tally verified, since kurtosis excess is quite ubiquitous. In 
this case, and since we are considering leptokurtic distri- 
butions, P (v) is also leptokurtic. An illustration of this 
sort of example is presented in Fig. [5] where a uniformly 

varies between and i. 

5 

Moreover, extending our range of values for exponent 
a, we can consider positive and negative values as we 
present in Fig. [5] For this last case, in the long-term, the 
system endures both platykurtic and leptokurtic regimes. 



For this case an analytical form is not, in principle, pos- 
sible to obtain. Nonetheless, we obtain the numerical 
solution for P (v) as we show in Fig. [7] 



B. A generalised Weibull distribution 

In this subsection we treat, the standard case where 
in Eq. ([8]) u> evolves on a superstatistical fashion rather 
than a. For this particular case, when fl = u>~ 2 follows 
a x 2 -distribution, 



1 



2+i 



exp 



v n 

"2 



(23) 



the long-term stationary distribution, P (v) , that corre- 
sponds to a weighted average of p (v) over all possible 
values of yields 



?w 0) = y< exp ? 



where, 



4 + i> — 6a — 2v a 
2 + v -2a-2va'' 



(v 2 )\ (24) 



v(\ -2a) fi 

(2 + i/-2a(l + i/)) 7 ' 
(25) 
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FIG. 5: Left panel: The full line represents the numerically obtained probability density function P («) vs. v for the case of 
superstatistical processes with uniform distribution, Eq. (|21l) . where ao — 0, and ai = § (7 = 100, a; = 10). The dashed line 
represents the PDF of the upper bound a = i and the dotted line the PDF of the lower bound a = (a Gaussian). The 
long-term PDF decays slower than the Gaussian. The grey symbols represent the PDF obtained from the numerical simulation 
of a superstatistical system with the same parameters 7 and to and PDF p (a). In the inset we show p (a) of that process. Right 
panels: Excerpt of a superstatistical time series of v which evolves according to Eq. (|8j) with 7 = 100, ui = 10, and a associated 
with a uniform distribution between a = and a = |. The time scale of updating a is 1 time unit which is rather larger than 
10~ 2 that is the time scale of relaxation towards stationarity. For this case «o = 0, and Qi = 4, with = 0.285 . . ., and 
kurtosis, k — 5.272 .... A fair similar distribution, namely Fig. 6(b), has been obtained in Ref. [2(j] for the velocities PDF of a 
long-range Hamiltonian system at a quasi-stationary state. 




FIG. 

P(v) 
form distribution, 



6: Numerically obtained probability density function 
vs. v for the case of superstatistical process with uni- 
Eq. (|21|l . where ao = — 5, and ai = g 
(7 = w = l). In this case v\j — 0.476..., and kurtosis, 
k = 3.129.... Interestingly, this example suggests the ex- 
istence of, at least, one distribution p(a) with ao / ai / 
for which P (v) is mesokurtic. 



and 



1 - 2 a. 



-2 a. 



(26) 



Owing to its comparability to the Weibull distribution 
inside S q framework, and following Ref. [2l| we shall call 
PDF (gl]) q- Weibull distribution. The q- Weibull distri- 
bution has been used to numerically adjust a variety of 
PDFs, but, to the best of my knowledge, no dynamical 



FIG 

P(v 

distribution 



7: Numerically obtained probability density function 
vs. v for the case of superstatistical processes with a x 2 - 
Eq. (|22p where v — 5, and a = 7 = ui = 1. For 



this example we have — 1.150 . . ., and kurtosis, k — 1.999. 



basis has been presented so far. Distribution (f2"4"|) belongs 
to the Burr class of probability density functions [22| . For 
small values of q W (v) in Eq. (j2"4")l goes to zero as a 
power law with exponent 6, and for large the same 
distribution also vanishes as a power law but with expo- 
nent a/ (1 — q) + b. For this case, it is possible to evaluate 
even moments of order m when the following conditions 
are verified, 



a > 0, b > 



b > 



(27) 
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When b — 0, i.e., a — 0, P (v) turns into a g-Gaussian 
distribution recovering the scenario of Ref. @]. In FigM 
we depict examples of q W (v) for some values of a, b, and 
<?■ 

IV. REMARKS AND PERSPECTIVES 

In this manuscript we have discussed a superstatisti- 
cal approach of a stochastic process which belongs to 
the Feller class, and whose (local) stationary probability 
density function is reminiscent of a Weibull distribution. 
Such proposal, which might be seen has a dynamical ana- 
logue for Ref. [Hj], has been driven on two ways. The 
first one in which the power of multiplicative noise term 
is time-dependent, and another one that is equivalent to 
the evolution of the noise width in the stochastic term of 
the corresponding differential equation. The main advan- 
tage of the former is the possibility of mimicking systems 
which hold a rather complex dynamics that goes through 
a (random) sequence of dynamical regimes. We have de- 
termined cither analytically or numerically the long-term 
probability density function which can assume all types 
of kurtosis. Concerning the latter superstatistical ap- 
proach, we have verified that it allows the emergence of 
a generalisation of the Weibull distribution within the 
framework of Tsallis non-additive entropy. This distri- 
bution has proved to be valid for numerical adjustments 
in a large variety of systems. For the q- Weibull, we have 
also been able to obtain platykurtic, mesokurtic, and lep- 
tokurtic distributions. In respect of is actual implemen- 
tation to the analysis of experimental data we must refer 
that it requires the development of techniques to capture 
both of p (a) and T as it was introduced for cases with 



fluctuations in u> [H 0| • The solution for this non-trivial 
challenge either on an experimental or theoretical level is 
certainly welcomed. 

Last of all, the extension of the idea of a stochastic 
dynamics with superstatistical exponents that lead to 
a consistent emergence of leptokurtic long-term proba- 
bility density functions by means of considering, e.g., a 
non-zero mean or the combination of independent multi- 
plicative and additive noises will be the focus of further 
research. 
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